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This work discusses gravothermal catastrophe in astro-physical systems and pro- 
vides an analytic collapse solution which exhibits many of the catastrophe properties. 
The system collapses into a trapped surface with outgoing energy radiated to a future 
boundary, and provides an example of catastrophic collapse. 
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I. INTRODUCTION 



■ n n 

Since 1968, Lynden-Bell [l 3 1 has illuminated the concept of negative heat capacity in 
astrophysical systems. He explained Antonov's theorem which states (roughly) that a 



> , 

spherical collection of self-gravitating point masses has no global entropy maximum. Ex- 
planation was needed since the most probable state for such a system is at the maximum 
qq , of the Boltzmann entropy. Lynden-Bell and Wood considered a self-gravitating gas sphere. 

They calculated the energy E and heat capacity Cy of the isothermal sphere. A graph of 
the sphere's binding energy vs density contrast showed a local maximum near an inflection 
point, and the plot of specific heat vs density contrast showed a stable branch for negative 
specific heat. 

Negative Cy can be understood as follows: the virial theorem for inverse square forces 
of bounded systems relates the average kinetic energy and average potential energy as 

< K.E. > = -- < P.E. > . 
2 

With total energy E = < K.E. > + < P.E. > one has 

E = — < K.E. > negative, 
but for moving particles K.E. = ^NksT. It follows that the specific heat is negative 

dE 3 *r, 

C v = — = — Nk B . 
v dT 2 



2 



A negative Cy system in contact with a large thermal reservoir will have fluctuations that add 
energy and make its transient temperature lower, causing inward heat flow which will drive 
it to even lower temperatures. Thus negative Cy systems cannot reach thermal equilibrium. 

Lynden-Bell coined the name " gravothermal catastrophe" for stellar systems undergoing 
such collapse. To describe the collapse, we quote j^j]: "Conductive transfer of heat from 
the central region will raise the high central temperature faster than it raises the lower 
temperature of the outer parts. No equilibrium is possible; the center continues to contract 
and get hotter, sending out heat to the outer parts." 

Astrophysical systems where gravothermal catastrophe may occur [g] are older globular 
clusters with compact cores, and bright elliptical galaxies with high central density profiles. 
Lynden-Bell and Eggleton studied the core collapse of globular clusters. In particular, they 
calculated self-similar collapse features of a gravitating gas sphere. One of the things they 
learned from their study was, that to model the formation of a central black hole in a globular 
cluster, they needed to go beyond self-similarity and study more dissipative collapse. 

In this work, we present an analytic solution [7] with shear and radial heat flow (SRH) 
which models many of the features of catastrophic collapse. The system collapses into a 
trapped surface with outgoing energy radiated to a future boundary. 

In the following section, the SRH collapse metric with pressure and density is developed. 
The results are given in section III and summarized in section IV. This is followed by two 
appendices. (A) contains the original SRH solution with arbitrary functions, and (B) has 
mass and trapped surface equations. 



II. COLLAPSING SRH FLUID 

I — I 

The SRH spacetime (7j is divided into an exterior region covered by the Vaidya metric 
and a collapsing interior region with spherically symmetric metric 

g^dx^dx" = A 2 dt 2 - B 2 dr 2 - R 2 dQ 2 (1) 

where A = A(t,r), B = B(t,r), R = R(t,r), and dQ 2 is the metric of the unit sphere. The 
stress-energy tensor is given by (G — c — 1) 



(2) 
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p is the isotropic pressure, w is the mass-energy density, ^ = g^—u^u", and q 11 is the radial 
heat flow vector orthogonal to u^. Time is comoving with ii^d^ = A~ l d t , u^u^ = 1. We 
use notation of Taub [8[ for the mass-energy density w. Taub's purpose was to distinguish 
mass-energy density from proper mass-density p, where w = p(l + e) and e is specific internal 
energy. 

The original SRH solution, given in Appendix A, provides the following metric functions 
and physical scalars (overdots denote d/dt and primes denote d/dr). 



R(t,r) = e t / t0 [(/3 2 /2)e- 4 ^ + (r + r ) 2 ] 1/2 , 

A(t, r) = aoR, 

B(t,r)=p /R. 



(3a) 
(3b) 
(3c) 



Metric functions A and B should be dimensionless and R should have dimensions of length. 
We write A = ao{R/c) so it is clear that «o is dimensionless. The other constants have the 
following dimensions: dim(to) = time, dim(/?o) = length, dim(ro) = length. 

The heat flow vector and scalar are given by kitq^dx^" = Qdr. The heat flow scalar is 



Q=--C-^y t/to . 

a R- 5 



The pressure is 



87cp 



(r + r ) 2 e 4t/t ° 



R 2 + /3 2 e- 2 */* 



(3 2 R 2 [R 2 - /3 2 e~ 2 */*°_ 



(1 + 1/ag) 
R 2 ' 



and the mass-energy density is 



1 

R 2 



Of 



(r + r ) 2 At/tQ 

Pi 



PI 



,2t/t 



w is negative when 



o; 



_ x _ &l{r + r ) 2 e4t/tQ 



2a R2 _ e 2t/t H < q 



a 2 (3 2 - (3 2 - a 2 (r + r ) 2 e 4 ^ - 2a 2 R 2 e 2t ^ < 
or ft* + 3a 2 (r + r ) 2 e 4t/to > 0. 



(4) 



(5) 



(6) 



This inequality always holds, and so w is negative. 
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III. COLLAPSE RESULTS 
A. End stage of collapse 

For small {3q, R goes as 



i?~e* /to (r + r ). (7) 



The pressure and mass-energy density, at late times and for small (3q, go as 

e 2t/t 

e 



8np ~ (8) 

™ - - 3r ^- ( 9 ) 

Po 



At late times, the magnitude of the mass-energy density increases exponentially, and the 
fluid approaches a photon gas with equation of state 

p= - | w | . (10) 

Negative mass-energy density w is linked to the gravothermal catastrophe. The collapsing 
fluid has a negative specific heat which derives from negative internal energy e. While 
both w and e are negative, the proper mass density, p = w/(l + e), remains positive. To 
separately compute e would require a complete thermodynamic analysis. One would need a 
causal relativistic description of thermodynamics such as the Israel- Stewart 'second order' 
type theory 0, 1(3]. This has been left for future work. 



As a model of gravitational dyamics, Chavanis et al [ll| studied the collapse of a gas 
of self-gravitating Brownian particles in a closed sphere. For catastrophic collapse in the 
microcanonical ensemble, they found the mass density to go as 

p ~ p (r /r) 2 - 21 (11) 

and further, in a stability study, they found the perturbation Sp/ p has a " core-halo" structure 
(hinting at central black hole formation). Without developing the thermodynamics of a 
collapsing SRH fluid, the density approximation in Eq. ( flTT) can be used to compute internal 
energy from e = w/p — 1. With w at late times given in Eq.flH]), we find 

e 2t/t 

e ~ -const — -j (r/r ) 2 - 21 - 1 (12) 
PoPo 

When the density is given a temperature profile (in the canonical ensemble p ~ T -1 / 2 ), then 
e = e(T) yields negative specific heat cy = de/dT. 



5 



B. Trapped Surface 

From Eq. (1B2j) the Misner-Sharp mass within a t = const, r = const 2-surface is 



2m = R 



1 _ (r + r ) 2 /tQ 
«0 Po 



(13) 



Null rays entering and leaving the 2-surface are described by the expansions of their respec- 
tive generators. Equations flB6|) and flBTj) show that both null generators have non-negative 
expansions and so a trapped surface exisits at R = 2m. At the trapped surface, Eq. (I13p sets 
a value for constants a and (3 . 



4 = Kap + r ) 2 e 4 ^/ i ». (14) 



The expression for R, Eq.f l3"a|) . yields 

< a p = (rt rap + ro)/3o^^. (15) 



We see from Eq. (HI) that, at late times, the heat flow scalar goes to zero 

Q — — — — : — ~TT777 y 0- (16) 

a (r + ro) 2 e'/'° 

The heat flow shuts off as the fluid collapses into the trapped surface. 



C. Rate of collapse 

The rate of collapse scalar is 6 = V^m^ = — l/(a R). We again quote Lynden-Bell: 
"During the gravothermal catastrophe ... the center continues to constrict and get hotter, 
giving out heat to the outer parts, but the temperature difference increases and drives the 
collapse onwards still faster." At distances just beyond the trapped surface 

a f3 Q ^ ^ 

At late times the rate of collapse increases exponentially. This is a necessary component for 
a model of catastrophic collapse. 

IV. SUMMARY 



An analytic solution of Einstein's equations for dissipative collapse has been presented. 
The original SRH solution contains arbitrary functions of time which have been chosen here 
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to provide an explicit solution. The system collapses into a trapped surface with outgoing 
energy radiated to a future asymptotic boundary. The collapsing fluid has negative mass- 
energy, which has been related to negative specific heat. Lynden-Bell has linked collapse with 
negative heat capacity to gravothermal catastrophe. The SRH collapse has many features 
that model the gravothermal catastrophe. 



Appendix A: Original SRH Solution 

This solution is given in equations (14), (15), and (16) of j3] [there is an error in the 
first term of Eq.(16). fthi should be (/3g/2)/ii]. The metric components are, with arbitrary 
functions hi(t) and h 2 (t) 



A(t, r) = a R, 
B(t,r)=p /R, 

R(t,r) = [(ft/2)h 1 + h^(r + h 2 ) 2 ]^ 2 . 

The match of g SRH to exterior Vaidya was done in [71 . 

The pressure, with parameters and arbitrary functions unchosen, is 



a R2 



, R' x 
-~R' 



+ 2- 



R^R_ 
R R 



1 



(R'f + (R 2 yS- 

R 



1 




a 2 R 2 









1 



(r ±M )2 



+ 



2(^)[(^)-a 



PlR 2 (3 2 [R 2 - ffih! - 2(h 2 /k)(r + h 2 )} 
with mass-energy density 

1 



(i + i/qg) 

R 2 



8lTW 



alR 2 



a 2 R 2 



"R } 



I? 

?2 



On 



ft 

ot 2 {r + h 2 ) 2 



2— + 3(-) 2 
R y R ! 



1 

R 2 



The heat flow vector and scalar are given by Anq^dx 11 = Qdr, 

a R? 

r + h 2 



Q 



a R 2 



(Ala) 
(Alb) 
(Ale) 



(A2) 



(A3) 



(A4) 
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For the particular case above, we choose 

hl = e -*/to 7 h 2 = r . (A5) 

Appendix B: MASS and TRAPPED SURFACE 

The collapse metric is spanned by the tetrad 

u^dx^ = Adt, r^dx^ = Bdr, 
h^dx^ = Rd'd, (p^dx 11 = Rsmtfdcp 

such that 

9iw = u^Uu - r il f v - - (Bl) 
The unique mass [3] m within t = const, r = const 2-surfaces is 

2m = R 3 R a ^J a ^^ip v = R{1 + R 2 /A 2 - {R'f/B 2 \. (B2) 

The metric is Petrov type D and the only non-zero Weyl tensor component ^2 is invariantly 
expressed as 48(\I/2) 2 — CaPni/C • The two principal null vectors of metric (GQ), normal to 
(■&, ip) 2-surfaces, are 

Pfy = a-% + B~ x d r (B3) 
■nPdp = A- l d t - B~ l d r (B4) 



with respective expansions 



R { A + B h R [ A B Y (B5j 



Consider a 2-surface S generated by Z M and n/*. If both null generators converge, i.e. ex- 
pansions 1? and have the same positive sign, then the 2-surface is trapped [13]. When 
R = 2m, Eq. (1B2|) provides 

R_R_ 
A~~B 

which implies 



•4=°. <r,=^- < B6 > 



Using A = a R, the expansion of / M is 



K = (B7) 



When ao > both expansions are non-negative signifying a trapped surface at i? = 2m. 
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